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The dynamical regimes of models belonging to the Kardar-Parisi-Zhang (KPZ) universality class 
are investigated in d = 2 -I- 1 by extensive simulations considering flat and curved geometries. 
Geometry-dependent universal distributions, different from their Tracy- Widom counterpart in one- 
dimension, were found. Distributions exhibit finite-time corrections hallmarked by a shift in the 
mean decaying as , where /3 is the growth exponent. Our results support a generalization 
of the ansatz h — Voot + {Tt)^x + ^7 + C,t^^ to higher dimensions, where Voo, T, C, and rj are non- 
universal quantities whereas /3 and x ai'e universal and the last one depends on the surface geometry. 
Generalized Gumbel distributions provide very good fits of the distributions in at least four orders 
of magnitude around the peak, which can be used for comparisons with experiments. Our numerical 
results call for analytical approaches and experimental realizations of KPZ class in two-dimensional 
systems. 
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Almost three decades after Kardar, Parisi and Zhang 
(KPZ) yj proposed their celebrated equation to describe 
the coarse-grained regime of evolving surfaces, a renewed 
burst of interest on it has been stood out due the exper- 
imental realization of its universality class in turbulent 
Hquid crystal setup [2] and the achievement of invaluable 
analytical solutions for distinct dynamical regimes and 
geometries in d = 1 -I- 1 j^. The KPZ equation reads as 

^ = ^v^/^ + ^(v/.)^ + C, (1) 

where ^ is a white noise of mean zero and amplitude \/lJ. 
Despite of its original conception for evolving interfaces, 
the KPZ equation has also found its place in others im- 
portant physical systems [1]. 

A great advance in the theoretical understanding of 
the KPZ universality class has begun at early 2000s with 
the seminal works of Johansson [5 and Prahofer and 
Spohn 6^ presenting analytical asymptotic solutions of 
some models in the KPZ class. These solutions link 
the height's stochastic fluctuations to universal distribu- 
tions [7J of the random matrix theory. Inspired in these 
exact results, the ansatz 

h = v^t + sx(Ttfx. (2) 

with the exactly known growth exponent jS = 1/3, 
was conjectured as describing the asymptotic interface 
fluctuations of any model belonging to KPZ class in 
(i=l + l|4l|8]. In this equation, s\ = sgn(A), while 
the asymptotic velocity Voo and T are model dependent 
parameters and x is a universal random variable with 
time-independent distribution given by the Gaussian or- 
thogonal ensemble (GOE) for flat geometries [3 E] and 
the Gaussian unitary ensemble (GUE) for the curved 
ones jSlll]. Notice that, in terms of the constants of the 
KPZ equation, the parameter T is given by F = |A^|A| 



with A = D/v jl]. These geometry-dependent univer- 
sality were confirmed in turbulent crystal liquid experi- 
ments [2 and in stochastic simulations of several models 
without known analytical solutions [5tillj. 

Many fine-tuning results have been aggregated to 
the asymptotic height distributions (HDs) of one- 
dimensional KPZ systems. The limiting processes de- 
scribing the surface fluctuations are known as Airyi 
and Airy2 processes for flat [T2 [T3] and curved geome- 
tries [14] , respectively. Finite-time corrections to Eq. ^ 
were also analytically [31 [TS] , experimentally [2j , and nu- 
merically [12] observed, leading to the generalization 

h = v^t + sx{TtYx + V + C,t'^. (3) 

where 77 and C are non-universal. The correction 77 in- 
troduces a shift in the distribution of the scaled height 
q = i'^"(rt)* relation to the asymptotic distributions. 
The hallmark of this correction, a shift in the mean van- 
ishing as [q) — (x) ^ t^^^^, has been verified in the crys- 
tal liquid experiments [5] and computer simulations of 
several models [HI UHl HI- To our knowledge, only two 
exceptions have been reported. In the first one, Ferrari 
and Frings [15] analyzed the partially asymmetric sim- 
ple exclusion process and found a specific value of the 
asymmetry parameter where there is no correction up 
order 0{t^^^^). Off-lattice simulations of an Eden model 
consistent with a decay t~^/^ were reported [TT], but a 
subsequent analysis showed that the unusual behavior is 
an artifact of low precision estimates of f 00 and a long 
crossover to the scaling law t~^/^ [l6] . 

In contrast to the deep understanding of the KPZ class 
in c? = 1 + 1, essentially no exact results are available 
in d = 2-1-1, the most important dimension for ap- 
plications |17j . Indeed, available analytical approxima- 
tions [H] fail in predicting the best numerical estimates 
of the scaling exponents ITHj. The best we know about 
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the KPZ class in d = 2 + 1 comes from simulations: 
The scaling exponents [19] and height distributions in 
the stationary regime [20] are accurately known and its 
universality has been verified. A few works, impaired by 
finite-size effects, had studied height distributions in the 
dynamical regime using flat geometry |2H 122] when, very 
recently, Halpin-Healy [SS] reported large-scale simula- 
tions of some KPZ models that convincingly suggest the 
universality of the height distributions. Halpin-Healy's 
analysis is in consonance with our results. 

In the present work, a detailed study of the dynam- 
ical regime of several KPZ models in 2 -I- 1 dimensions 
is presented. Both flat and curved geometries are con- 
sidered. We go beyond the Halpin-Healy's results and 
show that the generalized KPZ ansatz given by Eq. ([3| 
still holds in d = 2 + 1 with the proper growth expo- 
nent /3 = 0.24. The universality of x, that differs from 
the counterparts in H-1 dimensions, is confirmed and its 
geometry-dependence characterized. Also, we have veri- 
fied that the corrections in the mean vanish as and 
non-universal corrections were found for higher order cu- 
mulants. We compensate the absence of an exact ana- 
lytical expression for the HDs, showing that generalized 
Gumbel distributions [23] fit noticeably well the heights 
scaled accordingly Eq. ([s]). 

Flat geometry - We performed extensive simulations 
of three models in the KPZ class, namely, the restricted 
solid-on-solid (RSOS) single step (SS) [II] and etch- 
ing |26j models. Square lattices with up to 2^^ x 2^^ sites 
and periodic boundary conditions were used. Except for 
SS model, for which a checkerboard initial condition was 
used, an initially smooth substrate was considered. Up 
to 10"^ runs were used in averages. 




FIG. 1. (Color online) Main plot: Height variance against 
time for RSOS and SS models. Dashed line represents the 
power law Bottom inset: Effective growth exponents 

against time. The horizontal line represents the accepted KPZ 
value in d = 2 -|- 1. Top inset: Interface velocity against f''"^. 

The growth exponent can be determined from w"^ = 
{h?)c ~ t^^ , where {X")c represents the nth cumulant 



of X. Figure [T] shows the evolution of the variance for 
two models, while the corresponding effective growth ex- 
ponents (local slope in curves Inw vs. Int) are shown in 
the bottom inset. The growth exponents obtained for all 
models are shown in Tab. |TJ in which an excellent agree- 
ment with the accepted KPZ exponent f3 = 0.24 is ob- 
served for all flat models. Differentiating (h) in Eq. ([s]) 
one finds {h)t = Voo + s\PT^ {x)t^~^ -I- ... A hnear re- 
gression in {h)t against t^^^ for t — > cx) yields Vryo- This 
procedure is illustrated in the top inset of Fig. [T] and the 
estimates for all investigated models are given in Tab. [T] 
The quantity P*^ (x) can be obtained from the asymptotic 
value of gi = {{h)t — Voc,)t^~^ / (3. It was shown that the 
value of P determined from gi is more reliable than using 
cumulants of order n> 2 [TB] . 




FIG. 2. Top: Determination of mean shift {rj} for fiat and 
curved geometries. Dashed line is the decay Bottom: 
Normalized cumulants g* = gn{t) / g„{oo) against scaled time 
for SS (left) and RSOS (right) models in a flat geometry. Lines 
are (scaled) linear regressions used to determine (7„(oo). 

The accuracy in determining universality in simula- 
tions may be very sensitive to the correction rj de- 
pending on the model and the attainable simulation 
time. So, it is important to determine the strength of 
corrections before analyzing height distributions. The 
mean (77) can be determined using the height scaled 
in terms of directly measurable parameters Voo and gi 
as q' = {h — Voat)/{s\git^) [TB]. Equation ^ implies 
1 - {q') = -{sx{'n)/gi)t-'^ + ... Figure [2] shows that the 
power law t~^ describes very precisely the shift, analo- 
gously as observed ind=I + I[lllllini[ISl[I3- So, 
using the prefactor of the power law t~^ , we determined 
{rj) for all investigated models. The estimates are shown 
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model 




51 


52 


53 


54 


iv) 




R 


S 


K 


RSOS 


0.31270(1) 


-0.773(1) 0.1936(4) 0.0364(3) 


0.0130(5) 


-0.5(1) 


0.240(3) 


0.324(3) 


0.427(5) 


0.347(8) 


SS 


0.341368(3) 


-0.881(1) 


0.250(1) 


0.0536(3) 


0.0219(5) 


-0.4(1) 


0.239(5) 


0.322(2) 


0.428(5) 


0.35(1) 


Etching 


3.3340(1) 


-2.348(3) 


1.715(3) 


0.950(2) 


1.00(1) 


0.6(1) 


0.235(5) 


0.311(2) 


0.423(2) 


0.340(5) 


RSOSC 


0.3134(2) 


-2.116(2) 


0.272(2) 


0.0481(6) 


0.0158(5) 


-1.7(1) 


0.24(1) 


0.061(3) 


0.339(8) 


0.21(1) 


SSC 


0.12611(2) 


-0.797(1) 


0.051(2) 


0.0037(2) 


0.00053(8) 


-1.2(1) 


0.23(2) 


0.080(5) 


0.32(4) 


0.20(5) 


Eden (001) 


0.6495(3) 


-3.543(3) 


0.785(8) 


0.234(3) 


0.13(1) 


9.8(5) 


0.243(7) 


0.063(2) 


0.336(9) 


0.21(2) 


Eden(lll) 


0.6242(2) 


-3.219(5) 


0.610(8) 


0.164(3) 


0.083(5) 


8.8(5) 


0.239(6) 


0.059(2) 


0.34(1) 


0.22(2) 



TABLE I. Non-universal and universal quantities for the dynamical regime of KPZ models. Definitions in the text. 



in Tab. [H 

From Eq. (§, we have that scaled cumulants (?n(^) = 
(/i")c/(s"i"'^), n > 2, converge to T"'^{x")c for t oo. 
Contrasting with the first cumulant, the corrections in 
Qn depend on the model. Figure [2] bottom shows the 
scaled cumulants against t~^^ where A was assumed in- 
teger (used values are indicated nearby each curve). For 
sake of visibility, curves were normalized by the asymp- 
totic value r"'^(x")c obtained by extrapolation in plots 
gn versus t~^^ . These estimates are shown in Tab. |l] 
For SS (bottom left in Fig. [2| and Etching (data not 
shown) models, the corrections are quite consistent with 
(9")c — {x^)c ~ t""'', in analogy to the exact solution of 
the KPZ equation with edge initial condition and exper- 
imental results in d= 1 + 1 [2l[3]. However, in RSOS the 
second cumulant present a different behavior with the 
shift decaying approximately as t~^^ demonstrating the 
non-universality of the corrections in cumulants of order 
n > 2. 

The parameters gi, i = 1 to 4, shown in Tab. |l] depend 
on F, which can not be determined directly from height 
distributions [16 . However, one can investigate dimen- 
sionless cumulant ratios that are independent of F and, 
therefore, are expected to be universal. In Tab. ^ we 
show the ratios R = g-i/gl = (x^)c/(x)^ S = g^/gj"^ ^ 
{X')c/{x^fj^ (skewness) and K = g,/g^, = {x^)c/{x')l 
(kurtosis) for all investigated models. The ratios for dif- 
ferent flat models are essentially the same, confirming the 
universality of x conjectured initially. Notice that they 
are different from the ratios for GOE distributions ex- 
pected for their one-dimensional counterparts [Ij. Since 
an infinite hierarchy of cumulant ratios can be measured, 
in principle, we can determine all cumulants in terms of 
the first one. Our estimates for S and K are in good 
agreement with those found by Halpin-Healy in |23) , but 
fluctuating estimates for (x) and (x^)c presented there do 
not allow a reliable estimate of R (values ranging from 
0.33 to 0.51 are extracted from Ref. [23]). We believe 
that the corrections in distributions, mainly in the mean, 
are responsible by the apparent non-universality of R in 
Ref. |55]. Our estimates of S and K are also consistent 
with former, small-size simulations [21 and also with re- 



cent simulations of Eden model on flat substrates |27) . 
confirming the universality of the HDs. 

Due to the lack of rigorous results in 2-1-1 dimensions, 
we are currently not able to associate our numerical re- 
sults to an analogous of TW distributions. However, 
previous works dealing with linear systems have shown 
that the generalized Gumbel distribution with a non- 
integer parameter m fits the probability density func- 
tions of stationary quantities in several equilibrium and 
non-equilibrium systems [SH [28l [29] . We have obtained 
a very good agreement between our simulations and the 
so-called Gumbel's first asymptotic distribution of mean 
{X) and variance (X^)^ [^, 

77? ^ h 

G(X;to) = — -exp[-m(zx + e-^-)], (4) 
i (to) 

where b = VWMTC^^: = ^((-^> - X + s), 
s — [Inm — ij^o{m)]/b, F(X) is the gamma function and 
i>k{X) the polygamma function of order k [30]. The pa- 
rameter m allows to change simultaneously, but not inde- 
pendently, the skewness and kurtosis of the distribution. 
For m = 6, one obtains a skewness Sq = —0.4247 and 
kurtosis Kq — 0.3597 very close to the universal values 
for flat models shown in Tab. HI 

The height distribution scaled to a mean 1, ac- 
cordingly the non-universal parameters, becomes q* = 
{h — Vaot — {r]))/{sxgit^), leading to a variance {q*'^)c = 
R. In top panel of Fig. [3j the scaled heights for flat mod- 
els are compared with a Gumbel distribution for m — 6, 
mean 1 and variance R = 0.32. A remarkable collapse is 
observed around the peak for at least four decades. From 
an experimental perspective, it is extremely hard to mea- 
sure distribution extremes with an accuracy comparable 
to our simulations. Hence, the Gumbel approximation 
is a useful reference to check the KPZ universality class 
in 2-1-1 dimensions. Notice that in a linear scale, simu- 
lations are indistinguishable from the Gumbel distribu- 
tion in contrast with the TW distributions that not even 
barely fit the distribution's peak as can be seen in inset 
of Fig. [3] Interestingly, the rightmost tail of the scaled 
distributions is well fitted by the scaled GUE distribu- 
tion Xgue/ iXgue) ■ It is worth mentioning that general- 
ized Gumbel functions was compared with distributions 
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FIG. 3. (Color online) Height distributions for flat growth 
scaled to mean 1 compared with a Gumbel distribution for 
m = 6 and variance R = 0.32. The inset shows the same data 
in a linear scale. The growth times are: t = lO"* (RSOS), t = 
8000 (SS), and t = 2000 (etching). Scaled TW distributions 
are included for sake of comparison. 



FIG. 4. (Color online) Height distributions for curved growth 
scaled to mean 1 compared with a Gumbel distribution for 
m = 9.5 and variance R = 0.062. The growth times are : 
t = 1012.4 (RSOSC), t = 4000 (SSC), and t = 549.7 (Eden). 
Inset: Scaled height distributions disregarding the shift {rj}. 



of height extremes in the stationary regime of KPZ and 
other non-Unear models in Ref. [3T]. A good fit around 
the peak and large deviations in the tails were observed. 

Curved geometry - We study radial geometry using 
the on-lattice Eden D model [T6]. Due to the intrin- 
sic anisotropy of on-lattice Eden clusters, we investigate 
surface fluctuations along axial (100) and diagonal (111) 
directions. We also considered curved surfaces using the 
RSOS and SS models growing in a corner (RSOSC and 
SSC), where fluctuations in (111) direction are consid- 
ered. Details of the models and simulation are presented 
in Ref. [T5], where we carried out a detailed study in 
d = 1 + 1 and obtained the expected KPZ scaling, QUE 
TW, for curved growth. 

The growth exponents found for all models agree very 
weU with the KPZ value P = 0.24, as shown in Tab. |l] 
The non-universal parameters related to are shown in 
Tab. |l] The asymptotic velocity of SSC model has been 
under debate [321 133] and our estimate is in agreement 
with Ref. [33] ■ Again, the shift in the mean scales as 
exactly as in the flat case (Fig. [2]). However, the 
amplitude of the corrections are in general much larger 
than in the flat case, particularly for Eden model, and 
plays a central role for the time scale simulated in the 
present work. Corrections in gn are of order t~^^ or faster 
in analogy to the flat case. 

Dimensionless cumulant ratios are also universal for 
curved geometries as shown in Tab. [l] These ratios differ 
from those of the flat case and are even further from the 
TW values known for 1-1-1 dimensions. Our cumulant ra- 
tios are also in agreement with those reported by Halpin- 
Healy for a single model in the so-called point-point ge- 
ometry [33]. Once again, the scaled height distributions 



are well fitted by a generalized Gumbel distribution with 
m = 9.5, which has So = 0.335 and Kq = 0.224. A very 
important remark is that curves do not collapse if the 
correction (77) is not explicitly included in the analysis as 
shown in the inset of Fig. |4| Rescaling the distributions, 
accordingly to Eq. ([3]), to mean 1 and variance R = 0.062 
we found a good data collapse, with exception of the SS 
model (Fig. |4| . This is due to its larger value of R (pos- 
sible produced by large fluctuations). 

Assuming the last term in Eq. (|3| has the form Ct~^ , 
one has that 

s>.{{h)t - v^)t^-P = 31 - ^sx{C)t-^-^. (5) 

Our simulations show that gi converges to its asymptotic 
value with a correction quite close to t~'^^ in all fiat and 
curved growth models. So, the last term in Eq. ^ de- 
cays with an exponent j — (5. An equivalent result was 
obtained in the simulations of KPZ models in d = 1 -|- 1 
where a term t~^^^ was identified in the KPZ ansatz [16| . 
So, we have an additional evidence that the generalized 
KPZ ansatz in 0? = 1 + 1 has an equivalent counterpart 
in higher dimensions. 

In conclusion, we have studied the height distributions 
in the dynamical regime of KPZ systems in d = 2 -t- 1 and 
confirmed the universality of geometry-dependent distri- 
butions found very recently by Halpin-Healy [23]. How- 
ever, we have gone further and characterized also the 
finite-time behavior of the distributions. As in the 1-1-1 
case, the shift in the mean decays as t^^ and the cor- 
rections in higher order cumulants are non-universal and 
decay faster or equal than We also show that gen- 

eralized Gumbel distributions, commonly applied to fit 
distributions in linear systems [Ml |251 [53] , fit noticeably 
well the HDs of KPZ models that are non-linear. Such 
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distributions and the finite-time behaviors may play an 
import rule in the experimental study of KPZ systems. 
Furthermore, they may motivate and guide analytical in- 
sights to the understanding of the KPZ universality class 
in two dimensions. 

Authors acknowledge the support from CNPq and 
FAPEMIG (BraziUan agencies). 
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